Two-color nonlinear localized photonic modes 
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We analyze second-harmonic generation (SHG) at a thin effectively quadratic nonlinear interface 
between two linear optical media. We predict multistability of SHG for both plane and localized 
waves, and also describe two-color localized photonic modes composed of a fundamental wave and 
its second harmonic coupled together by parametric interaction at the interface. 
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Cascaded nonlinearities of noncentrosymmetric opti- 
cal materials have become an active topic of research 
over the last years due to their potential applications in 
all-optical switching devices Parametric interaction 
is known also to support solitary waves, in particular 
spatial quadratic solitons, which are two-frequency self- 
trapped beams consisting of a fundamental wave para- 
metrically coupled to its second-harmonic ||. Usually, 
solitary waves are considered for homogeneous media 
where spatial localization is induced by self-focusing and 
self-trapping effects. However, localized modes can exist 
even in a linear medium at defects or interfaces, and they 
are known as linear defect or interface modes. The prop- 
erties of nonlinear defect modes are usually analyzed for 
nonresonant Kerr- type nonlinearities ||. Here we con- 
sider a qualitatively different situation and introduce an- 
other type of nonlinear defect mode: a two-frequency ( or 
two-color) localized photonic mode, where the energy is 
localized due to the parametric wave mixing induced by 
an interface between two linear optical media or a thin 
layer with a quadratic (or x ) nonlinearity embedded 
in a linear bulk medium. 

The physical motivation for our model is twofold. First 
of all we point out a fundamental property of inhomoge- 
neous nonlinear optical media. Let us consider an inter- 
face between two semi-infinite bulk optical media, which 
are either clamped together or separated by an infinitely 
thin layer. If the bulk medium has inversion symme- 
try, then its quadratic nonlinearity must vanish. How- 
ever, the interface breaks the symmetry and therefore 
the interface nonlinearity should possess a nonvanishing 
quadratic response, due to a nonzero contribution from 
the spatial derivatives of the electric field Q . 

Second, there exists a strong experimental evidence 
of SHG in localized photonic modes. For example, re- 
cent experimental results reported SHG in a truncated 
one-dimensional periodic photonic band-gap structure, in 
which a nonlinear defect layer was embedded. An en- 
hancement of the parametric interaction in the vicinity 
of the defect was observed, suggesting that SHG occurs 
in local modes, while being completely suppressed for 
other propagating modes. If the band gap of the peri- 
odic structure is wide, we can describe this SHG process 
by a model with a local quadratic nonlinear defect. 



The main purpose of this paper is to introduce and 
study an analytically solvable model for SHG in non- 
linear localized modes, based on the (induced or inher- 
ent) quadratic nonlinearity of an interface separating two 
(generally different) linear bulk media or a thin nonlin- 
ear defect layer embedded in a bulk medium. In particu- 
lar, we describe two-color localized nonlinear defect modes 
where the energy localization occurs due to parametric 
coupling at the interface. 




FIG. 1. Scattering of a plane FF wave (dashed) on a 
X^'-interface between two linear media 1 and 2. The gen- 
erated SH (solid) can be either propagating or localized. 

We consider a fundamental frequency (FF) wave prop- 
agating along the Z direction in a linear slab waveguide, 
as shown in Fig. [j] We assume that an interface (or 
defect layer) at x=0 possesses a quadratic nonlinear re- 
sponse, so that the FF wave can parametrically couple 
to its second harmonic (SH) at the interface Q. The 
coupled-mode equations for the complex envelope func- 
tions Ej(x, Z) (j=l,2) can then be written in the form 
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n 2 {x)E 2 +T 2 (x)E 2 = 0, 



where Dj are diffraction coefficients (Dj > 0). For the 
geometry shown in Fig. [I] and the approximation of an 
infinitely thin interface layer (valid when the width of the 
layer is much smaller than the FF wavelength), we take 
rij(x) = noj(x) + KjS(x) and Fj(x) = jjS(x), where jj 
are the nonlinearity coefficients, Kj account for the phase 
velocity differences in the layer and bulk materials, and 
noj(x) — Uqj, for x > 0, and noj(x) = n^-, for x < 0. 

In order to reduce the number of parameters we 
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normalize Eqs. (0) as follows: E\{Z) = u(z)/ ^7172, 
J5 2 (^) = v(z)/ji, o = D 2 /Di, oti = K1/-D1, and 
«2 = K2/D1, where z = Z/D\ is measured in units of 
D\. Then the coupled equations take the form 



Qu 11 

i-g- + ^2 +vi(x)u + 6 (x) (a\u + u*v) = 0, 
i- — h + V2 (x) v + S (x) (012V + u 2 ) = 0, 
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where Vj{x) 



riQj/Di, for x > 0, and Vj{x) 



n j/Di, for x < 0. If the mismatch 2n 01 
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is small then er=l/2 is a good approximation, which we 
use below in the numerics. The system (^|) conserves the 
Hamiltonian 
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5{x) [ai|u| 2 + y H 2 + Re(uV)] } dx, 



and the total power P = (|u| 2 + \v\ 2 ) dx for spa- 

tially localized or periodic solutions. 

Scattering problem. To analyze the scattering process 
we use that Eqs. (||) are linear for 1 / 0, and write the 
total field as a superposition of plane waves, 



u{x, z) 
v(x, z) 



aie -l\ 1 z+iq 1 x + bie ~t\ 1 z-iq 1 x. x < Q 



Cie -i\iz+iq+x. 

b2e -i\ 2 z-iq-x. x < Q 
C2e -^2Z+tq+x. X> Q^ 



X > 0, 



where ai, &i, and c\ are the amplitudes of the incident, 
reflected, and transmitted FF waves, respectively. Corre- 
spondingly, 62 and C2 are the amplitudes of the generated 
SH waves on both sides of the interface. The dispersion 
relations are then given by 
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(3) 



and the continuity condition at x—0 yields the relations 
a\ + b\ = ci and &2 = C2. Next, integrating Eqs. (||) 
over an infinitely small segment around the interface, we 
obtain the relations between the field derivatives on op- 
posite sides of the layer, 



-[du/dx] 



-a[dv/dx]_ Q 



(u*v + aiu)\ x=0 , 
(u 2 + a 2 v) | x=0 . 



(4) 



This gives the phase- matching condition 2Ai = A2, which 
is a general requirement for stationary propagation of FF 
and SH without energy exchange, and two algebraic re- 
lations for the amplitudes, 



-i (<7! + qt) ci + 2iq l 01 = c*c 2 + ajci, 
-ia [fl2 + it) C2 = c 2 + a 2 C2- 



(5) 



All unknown parameters qf, q 2 , Xj, bj, and Cj can now 
be expressed in terms of the amplitude a\ and transverse 
wave number q^ of the incident FF wave. 

The SH amplitude at the interface is determined from 
Eq. (§), 
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-c\j [a 2 + ia (q 2 + q%)] 
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The wave numbers are found from the phase-matching 
condition 2Ai = A2 and the dispersion relations ([3|), 
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These values can be either real or imaginary, correspond- 
ing to plane waves and waves that are spatially localized 
at the interface, respectively. Note that the sign of the 
wave numbers is fixed according to the predefined geom- 
etry of the problem (see Fig. [I]) . 

Multistability. Substituting Eq. (||) into Eq. (||) we ob- 
tain the characteristic equation for the FF wave intensity 
I ci 1 2 at the interface, 
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where ot\ — ai+i (<?f + qf) and 5 2 = ct2 +icr (q 2 + q£) ■ 
Equation (@) is cubic in |ci| 2 , and thus three different 
roots may exist for a given input intensity |ai| 2 , corre- 
sponding to three different values of the amplitudes at 
the interface, as shown in Fig. ||(a). This describes a 
multistable SHG process. 

To study the stability of these solutions we investigate 
the corresponding linearized problem, similar to the anal- 
ysis of a different problem in Ref. Q . It is convenient to 
chose the perturbation functions with profiles which re- 
main self-similar upon propagation in the inhomogeneous 
medium. Then it can be demonstrated that the growth 
rate for nonoscillatory instability modes vanishes at the 
turning points <9|ai| 2 /<9|ci| 2 = 0. We have verified nu- 
merically that the growth rate is positive on the branch 
with negative slope [dashed line, Fig. |](a)], meaning that 
the corresponding solutions are unstable. However, a rig- 
orous stability analysis of all nonlinear modes is beyond 
the scope of this paper. 

Let us consider the simplest possible example, in or- 
der to illustrate the characteristic physical properties of 
the system. We choose the case when the linear media 1 
and 2 on each side of the interface (see Fig. [I]) are iden- 
tical, i.e., vf = Vj. It then immediately follows from 

Eqs. (^) that qj 1 = qj. This means that the SH can exist 
in two different states, propagating or localized, whereas 
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the FF waves are always propagating. The localized SH 
state can only be observed if j3 > 0, and then onl y for 
FF wave numbers less than a critical value: q\ < \J (3/2, 
where (3 = 1v\ — i>i- Note that localization does not 
depend on the wave amplitudes. 

For fixed material parameters ctj and (3 the multistable 
SHG regime can only be observed in certain regions of the 
input parameters q\ and |ai| 2 , as illustrated in Fig. ||(b). 
Importantly, multistability can be found for both propa- 
gating (region I) and localized (region II) SH waves. For 
other values of the material parameters the multistable 
scattering might be observed for a single type of the SH 
waves, or only a single-state field configuration may be 
possible. 
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FIG. 2. (a) Dependence of incident FF intensity ai| 2 on 
FF intensity at the layer |ci| 2 , illustrating existence of mul- 
tistability. (b) Regions of multistability in (gi,|ai| )-space 
for ai=l, a2=0.15, and /3=0.02. The three SH solutions are 
propagating (region I) or localized (region II). Only one solu- 
tion exists outside regions I and II. The dashed line gi=0.09 
corresponds to the case shown in (a). 

In Fig. H we summarize the different SHG regimes in 
terms of the material parameters. The diagrams are pre- 
sented for a\ = ±1, but they are invariant to the scaling 
(3 — > (3/a\, a-i — > a2/|ai|, meaning that they can char- 
acterize the multistability regions for any a.\. 



generated SH waves can be either propagating or local- 
ized. However, the situation when all the FF and SH 
waves are localized is also possible. These two-frequency 
nonlinear localized modes are of significant physical in- 
terest. To find stationary solutions for these modes we 
assume the following conditions in the general scattering 
problem: (i) no incident plane wave, i.e., ai=0, and (ii) 
all transverse wave numbers are imaginary, — i/Jtj, 

where ^ are real and positive (as the wave amplitudes 
should vanish at infinity). Then the amplitudes at the 
interface are |ci| 2 = y^i + fi^ — a%\ la + ^ ) — a^l 
and C2 = fJ-i + (J>i — ct\. Note that here only one wave 
number is arbitrary, all others are determined by Eq. (Q). 
Such localized states can exist for any combination of 
material parameters. Some examples are presented in 
Figs. |(a) and |(b). 



For the symmetric case, where v- 



and fi 
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the total power P — \c\\ //ii + |ca| I Hi can be written 
as a function of fii on ly The dependence P(pi) has al- 
ways, for any values of the material parameters, a branch 
with positive slope. Under certain conditions a second 
branch with negative slope may appear. This branch 
corresponds to smaller values of fii aim larger values of 
the Hamiltonian H, as shown in Figs. |J(c) and [|(d) with 
dashed lines. Thus we may conclude that for two bistable 
states the one with lower fj,i (he., higher H and negative 
slope dP/dfii < 0) is unstable. For other values of the 
material parameters the power ranges corresponding to 
branches with positive and negative slope may not fully 
overlap, or there can even be a gap. 
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FIG. 3. Regions I-III in (/3, a2)-space for ai = ±1 where 
multistability can be observed. In region I all three SH states 
are propagating, in II they can be either propagating or lo- 
calized, and in III they are all localized. The marked point 
corresponds to the case presented in Fig. ^| 

Nonlinear localized modes. As mentioned above, in the 
scattering of a plane FF wave the transmitted FF and 
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FIG. 4. (a,b) Intensity profiles of two different localized 
states with the same power P=20 (FF - dashed, SH - solid), 
(c) Power P(fii) with two branches, (d) Hamiltonian H vs. 
power P. The unstable branch is depicted with a dashed 
line. Points a and b mark the modes shown in (a) and (b), 
respectively. For all the plots on = 012 = — 1, P = 2. 
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Let us consider the generation of a stable two-color lo- 
calized mode by launching a localized FF wave at the 
interface. As the first step in the analysis of the dynami- 
cal problem we consider a simplified case, assuming that 
both the amplitude and phase velocity of the FF pump 
wave at the interface are constant. This is a reasonable 
approximation, if the initial FF wave is close to a station- 
ary mode localized due to a linear phase detuning at the 
layer (characterized by ax), and the generated SH wave 
remains small (essentially an undepleted pump approxi- 
mation) . For such a case, the original system (|^) can be 
reduced to a single equation for the SH wave, 



. dv 



d 2 v 



-2i\iz\ 



dx 2 



0, (9) 



with the initial condition w(0)=0. In Eq. (0) the initial 
intensity of the FF wave at the interface \cf \ and its 
propagation constant Ai are fixed. 

In the case a2—0, an exact solution of Eq. (^J) can be 
presented in the form, 



v(x, z) 
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i(2A!+i/ 2 )C p ix 2 /4<TC 



The expression under the integral describes decaying os- 
cillations with the increase of z. Thus, the amplitude of 
the generated SH exhibits oscillations as the solution ap- 
proaches asymptotically a stationary two-color localized 
state. We have performed a number of numerical sim- 
ulations (using a fully implicit finite-difference method) 
with Gaussian initial profiles of the FF wave and found 
that in the general case the formation of localized modes 
is accompanied by the same kind of transitional oscil- 
lations, as shown in Fig. ||. Furthermore, we also ob- 
served switching from a perturbed unstable two-color lo- 
calized mode, such as that shown in Fig. |](a), to a stable 
one. In the evolution process some energy is radiated, 
and the power corresponding to the localized mode is de- 
creased accordingly. Thus, a stable localized mode can 
only be generated, if initial power is above the threshold 
Pth = min Ml P(^i). 




FIG. 5. Generation of a two-color localized mode (material 
parameters are the same as in Fig. ^). Left: intensities of the 
FF (dashed) and SH (solid) modes at the interface. Right: 
evolution of the excited SH. The input power is P ~ 28 and 
the final state is close to that shown in Fig. ^b). 



In conclusion, we have introduced an analytically solv- 
able model for SHG in localized modes and predicted 
the existence of two-color nonlinear localized photonic 
modes supported by parametric interaction at an inter- 
face. Some of the properties of two-color localized modes, 
such as stability, generation, and switching, show a re- 
markable similarity with parametric solitons in homoge- 
neous optical x' 2 ' media Q and their interaction with 
localized perturbations of the mismatch parameter [Q. 
We believe that these results open a new class of prob- 
lems in the theory of nonlinear wave propagation in inho- 
mogeneous media with resonant nonlinearities, and they 
should be useful to understanding the fundamental differ- 
ence between the effects produced by nonresonant Kerr- 
type nonlinearities and those induced by parametric wave 
interaction. For example, the analysis of localized modes 
due to a single defect is the first step towards the 
theory of nonlinear modes and gap solitons in nonlin- 
ear photonic crystals ^ and the dynamics of the defect 
modes in such materials. 
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